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Abstract 

While LDA+U method is well established for strongly correlated materials with well localized 
orbitals, its application to weakly correlated metals is questionable. By extending the LDA Stoner 
approach onto LDA+U, we show that LDA+U enhances the Stoner factor, while reducing the 
density of states. Arguably the most important correlation effects in metals, fluctuation-induced 
mass renormalization and suppression of the Stoner factor, are missing from LDA+U. On the other 
hand, for moderately correlated metals LDA+U may be useful. With this in mind, we derive a new 
version of LDA+U that is consistent with the Hohenberg-Kohn theorem and can be formulated 
as a constrained density functional theory. We illustrate all of the above on concrete examples, 
including the controversial case of magnetism in FeAl. 

PACS numbers: 71.15.-m, 71.15.Mb, 71.20.Be, 71.20.Eh, 75.10.Lp 



*Permanent address: Dept. of Physics, South Dakota School of Mines and Technology, Rapid City, SD 57701 



One of the most influential, from practical point of view, developments in the Den- 
sity Functional Theory (DFT) in the last two decades was the LDA+U method (see, e.g., 
Ref. [I]). This method includes the orbital dependence of the self-energy operators, missing 
from the Kohn-Sham potential, in a relatively crude, pseudo-atomic way, neglecting the fine 
details of the spatial variations of the Coulomb potential. On the contrary, the standard Lo- 
cal Density Approximation, LDA accounts for the spatial variation of the Hartree potential 
exceedingly well, but neglects the orbital dependence of the Coulomb interaction. 

There is one inherent ambiguity in the LDA+U method: In LDA, all electron-electron 
interactions have already been taken into account in a mean field way. The Hubbard Hamil- 
tonian that represents the underlying physics of the LDA+U method also incorporates a 
large part of the total Coulomb energy of the system. Simple combination of the LDA and 
Hubbard Hamiltonian thus leads to a double counting (DC) of the Coulomb energy, so one 
may want to identify those parts of the DFT expression for the total energy that correspond 
to the interaction included in the Hubbard Hamiltonian and subtract them. However, since 
the DFT Hamiltonian is written in terms of the total density, and the Hubbard Hamilto- 
nian in the orbital representation, one cannot build a direct link between the two. Second, 
even if it were possible, that would be undesirable. Spatial variation of the Hartree and 
the exchange-correlation potentials is very important. It would be unreasonable to subtract 
that out just because it has been already taken into account elsewhere in a primitive way 
(roughly speaking, UN 2 /2). Rather, one wants to identify the mean-field part of the Hub- 
bard Hamiltonian, and subtract that, leaving only a correction to the LDA-type mean field 
solution. 

This is not a uniquely defined procedure. Several recipes exist, and it has been appreciated 
lately || that the results of LDA+U calculations may depend crucially on the choice of the 
DC recipe. It should be noticed that while for strongly correlated systems the LDA+U 
ideology is at least practically established, in a relatively new area of applying LDA+U to 
moderately-correlated, metallic systems 0, |3], the situation is very far from clear. 

In this Letter we analyze the effect of different DC prescription on the LDA+U results in 
correlated metals. We also present a systematic approach to the DC problem, of which the 
existing recipes are special cases. Finally, we discuss which problems associated with this 
class of materials can, in principle, be solved within LDA+U, and which cannot. 

We use for our analysis the spherically averaged form of the rotationally-invariant LDA+U 
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due to Dudarev et al. 
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where {7 and J are spherically averaged Hubbard repulsion and intraatomic exchange for 
electrons with the given angular momentum I, n ma is the occupation number of the m-th 
orbital, o = ±1 is the spin index, and the superscript means that the double counting 
terms have not been subtracted yet. Here p a mm i is the orbital occupation matrix, N a = Tr(p CT ) 
and N = J2 a N a . 

To subtract from Eq. (JTJ) the DC term, one naturally starts with the first two terms in 
Eq. ([]]), i.e., the Hartree and the Stoner energies. Both are explicit functionals of the spin 
density, and are likely to be better described by LDA. To identify the DC part of the last 
term of the Eq. ([!]), which explicitely depends on n mcr , is less trivial; one needs to work out 
a "mean field" approximation to this term, that is, substitute Tr(p a ■ p a ) by some quantity 
x a that depends solely on total spin density. Czyzyk and Sawatzky [0] suggested that x a 
should be equal to Tr(p CT • p a ) in the limit of the uniform occupancy, p <T r ^i A = 8 mm >n a , 
and, consequently, x a = {21 + l)n„, where n a = N a /(2l + 1). This leads to the following 
corrections to the total energy and the effective potential: 

cr 
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Here AMF stands for "Around Mean Field" [0] and <5p^ m , = p^ m , - n a 5 mm '- 

For strongly correlated systems the limit of the uniform occupancy is not correct (in fact, 
it is not correct even in weakly correlated systems, due to the crystal field). Thus, it is not 
surprising that for the systems with strongly localized electrons the AMF functional leads to 
rather unrealistic results. This observation led []7|, |8| to another prescription, x a = (2Z + l)n (J , 
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which produces the correct behavior in the fully localized limit (FLL) where n ma = or 1. 



Most of the modern LDA+U calculations utilize one of these two functionals, although in 
real materials the occupation numbers lie between these two limits. 

In the AMF the LDA+U correction to the electronic potential, Eq. |2|, averaged over 
all occupied states, is zero. This is a possible way to define a mean field (cf. the Slater 
approximation to the Fock potential), but not the way used in the DFT. The latter is a 
mean field theory that produces the correct total energy, not the correct average potential. 
AMF and FLL represent the "DFT" mean field if all occupation numbers are all the same, 
or are all or 1, respectively. It is easy to show that (21 + l)n^ < Tr(p CT ■ p a ) < (21 + l)n a , so 
that AMF always gives a negative, and FLL a positive correction to the total energy, while 
the right (in the DFT sense) recipe should give zero correction to the total energy. That 
can be achieved by using a linear interpolation between the two extremes corresponding to 
AMF and FLL, x a = (21 + 1) (an a + (1 - a)n 2 a ), where < a < 1, and 

&EEEI +u =-^~ £ [TW- 5p°)-(2l + l)an a (l - n a )} 
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In the spirit of the DFT, AE^ +u = 0, so 

E CT Tr(^-^) 
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We emphasize that a is not adjustable, nor is it a formal functional of the charge density, 
but it is a material-dependent parameter (like, say, U itself), defined by the self- consistent 
occupation matrix. However, in practical calculations it is better to recompute a after 
each iteration, as the current value of p a mm , changes. Note that the total energy is given 
by the regular LDA expression that only implicitly depends on U and J via the changing 
density distribution; it is variational with respect to the charge density at a fixed a, but not 
variational with respect to a itself. The fact that this prescription is derived according to the 
DFT ideology allows one to formulate the proposed LDA+U functional (unlike the existing 
LDA+U functionals) as a constrained DFT theory || at a given a, with the constraint given 
by Eq. (||). (U — J)/2 appears then as a Lagrange multiplier. 

We have tested the proposed functional on NiO, a prototypical compound for LDA+U 
calculations (see, e.g., ||10|| ). Fig. [I] shows the band gap and the mangetic moment of NiO as 
a function of U at J = 0.95 eV for three different functionals (Eqs. 0, |^, and calculated 
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within the linear-muffin-tin orbital (LMTO) method in the atomic sphere approximation 
(ASA). The parameter a ~ 0.5 is almost independent on U. Accordingly, the results of our 
calculations based on Eq. |] for both band gaps and magnetic moments lie right between 
those for AMF and FLL calculations, and the effect of U is reduced compared to the FLL 



calculations. This is in accord with a known observation HTDfl that in NiO the FLL LDA+U 



gives the best agreement with the experiment for U <6 eV, smaller than U ~ 8 eV calculated 
from the first principles] 
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FIG. 1: Mott-Hubbard band gaps, and magnetic moments of antiferromagnetic NiO for three 
flafors of LDA+U. The upper and lower values of the "error bars" correspond to the FLL and 
AMF functionals, respectively. 



Our next example is a weakly correlated metal FeAl. This paramagnetic material has at- 
tracted attention due to a recent suggestion by Mohn et al. 0] that the short-range Coulomb 
correlations within the LDA+U may be responsible for suppression of ferromagnetism found 
in all LDA calculations. More specifically, they found in their AMF LDA+U calculations 
a reduction of the density of states (DOS) at the Fermi level, Dp, which was sufficient to 
make the Stoner criterion smaller than 1 and stabilize the paramagnetic state. To analyze 
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this result, it is important to revisit the Stoner theory for the LDA+U case. 

In DFT, the Stoner parameter / is defined as I = —2d 2 E xc /dM 2 , the second derivative 
of the exchage-correlation energy with respect to the total magnetic moment. The param- 
agnetic ground state is unstable when Dpi > 1. This can be derived from the force theorem, 
which states that the total energy for small magnetizations can be computed by assuming 
a rigid shift of the bands by b = ±M/2Dp, so that the gain in the interaction energy, 
—IM 2 /A, competes with the loss in the one-electron energy, M 2 /4Dp. In the LDA+U the 
criterion holds, but the product Dpi changes, not only because Dp changes, but also be- 
cause the newly added interaction energy depends on M. Indeed, the force theorem calls 
for a change £p^ m / = baD mm i, where D mm > = — 7r _1 Im G mm /(Ep). When applied to the 
functionals Eq. (|2]) - Eq. (f|), it generates a change in the interaction energy, which results 
in an additional contribution to the Stoner parameter, 



(U—J) I (2Z+1). Given that the LDA Stoner parameter, J, is of the same order as J, we obtain 
for the total Stoner parameter I F ll ~ (U + 2U)/(2l + 1), which the well known expression 
for the Stoner factor in the atomic Hubbard model. On the contrary, AIamf (ot = 0) in 
this limit is zero. In real metals D mm i is complicated due to crystal field effects. Let us 
consider, for illustration, d-electrons in a cubic environment, and introduce the difference 
AD = D eg — D t 2 9 , where D eg and D t 2 g are e g and t2 g DOS per orbital at Ep, as a measure 
of the crystal field. This gives rise to a contribution to AIamp = ^{U — J)(AD / Dp) 2 . 
However, when LDA+U reduces Dp, and AIamf is not large enough to overcome the 
decrease in Dp, LDA+U may stabilize the paramagnetic state(cf. Ref. [0)) as ; f° r instance, 
observed in a very narrow range of large f/'s for FeAl by Mohn et al. p[ (of course, only in 
the AMF functional; the FLL functional produces a large AI » (U—J)/5, always increasing 
the tendency to magnetism). 

With this in mind, we performed LMTO-ASA calculations for all three LDA+U func- 
tionals, using fixed J = 0.95 eV. The results for [/-dependence of the magnetic moment and 
a are shown in Fig. |2| and compared with those by Mohn et al. 0. In our AMF calculations 
we also found a paramagnetic solution for U = 4.85 eV, which however coexsits with a 
ferromagnetic high spin solution (Fig. 0). Note that for well localized orbitals there is no 





In the limit of the uniform occupancy, Eq. (Q) for the FLL case (a = 1) reduces to 
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difference whether the (U — J) term is applied inside the atomic sphere or only inside the 
MT sphere, as in Ref. ||; however, in less localized cases, where a noticeable part of the 
o?-orbitals spills out of the MT sphere, the effect of the same U is smaller when applied only 
inside the MT sphere. One can see in Fig. ^that, indeed, our calculations with large U yield 
large a's and agree very well with Ref. |2J, while for small U (small a) the effect of U in our 
ASA calculations is stronger than in Ref. H. 
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FIG. 2: Magnetic moments of FeAl for AMF and DFT flavors of LDA+U compared with the 
results of Mohn et al. (Ref. §) 



All LDA+U functionals shift unoccupied bands up and occupied bands down. Therefore 
LDA+U broadens the bands crossing the Fermi level. Because of this broadening, in FeAl 
for small U the parameter a is initially decreasing (Fig. |2|) with a minimum a = 0.05 
at U = 2 eV. The magnetic moment also decreases in this region. At larger U, a starts 
growing again. At this point it is instructive to apply the logic of the constrained LDA 
approach in which for every fixed a the total LDA energy is minimized under the constraint 
^2 cr Tr(Sp a ■ 8p a )/[{2l + 1) J2a n v(^ ~ n o)\ = a i (U ~ -0/2 being the Lagrange multiplier. 
For a < 0.087 (Fig. |^) of the two possible solutions with U < 2 eV and U > 2 eV we should 
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choose the one with lower energy (smaller U). As a result, we find two admissible domains 
for U : an AMF-like with U < 2 eV and a FLL-like with U > 5 eV. The latter is clearly 
unphysical. Both solutions are ferromagnetic. The solutions with intermediate values of U 
and reduced magnetic moments are inadmissible in the framework of the constraint DFT 
formulation. 

On the contrary, our explanation of the paramagnetism in FeAl is that the ferromag- 
netism instability is suppressed by the critical spin fluctuations. There are many other 
systems for which the fluctuations in the vicinity of a quantum critical point reduce the 
tendency to magnetism. Further examples include Sr 3 Ru207 (Mlda ~ 0.8 /is, M exp = 0), 
ZrZn 2 [Mlda ~ 0.7 /ig, M exp = 0.2 /ig), and other. The physics that is missing from both 
LDA and LDA+U equations in such systems can be described as exchange of virtual elec- 
tronic excitations, roughly speaking, plasmons or (para)magnons. This leads to "dressing" 
of the one-particle excitations in the same way as the electron-phonon coupling "dresses" 
electrons near the Fermi surface, although in a correlated metal such mass renormalization 
effects occur on a large energy scale (of the order of U or J). LDA calculations cannot re- 
produce such a dressing, which has been observed in many different ways experimentally. 
For instance, LDA calculations do not explain large mass renormalizations in Sr 2 Ru04 flT3f . 
and large specific heat renormalization in many correlated metals, produce too large plasma 
frequencies, e.g., in YBa 2 Cu 3 7 , yield an optical absorption spectrum in Cr0 2 shifted by 
about 20% to higher frequency, as compared with experiment [TJ], and overestimate the ex- 
change splitting in Ni by a factor of 2 |15| . In all these cases the total width of the <i-bands 
is decreased, as opposed to broadening inherent to LDA+U. Here the essential physics is 
missing from the LDA+U as well as in LDA, while the spatial variation of the mean-field 
Coulomb interaction is treated better by the LDA. The missing physics is associated, to a 
large degree, with dynamic fluctuations. 

The dynamic version of the LDA+U method, the Dynamic Mean Field Theory (DMFT) 
16| , which can account for some spin fluctuations |]17l , resolves many of these problems. 



For instance, the mass renormalization in Sr 2 Ru04 is 3-4 fL3[ , far greater than possible 
renormalization due to the phonons. We applied all three flavor of LDA+U to Sr 2 Ru04 
and found no mass renormalization compared to LDA. On the other hand, Eliashberg-type 
calculations JTBJ of the renormalization due to spin fluctuations, using a spectrum deduced 
from the LDA band structure, give mass renormalizations of the right order. Similarly, 
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FIG. 3: FeAl density of states, D(E), in DMFT (solid line) compared with the nonmagnetic LDA. 
The DMFT solution is stable, the LDA is not (a ferromagnetic solution is stable), despite the same 
D(E F ). 



DMFT explicitely narrows the bands in Sr 2 Ru04 and enhances the electronic mass [19 



With this in mind, we applied the DMFT with a realistic U = 2 eV to FeAl and found 
the paramagnetic state to be perfectly stable, whith the bands narrower than in LDA, and 
the density of states practically the same (Fig. |). In other words, the spin fluctuations 
effectively reduce the Stoner factor /. 

To conclude, we observe that no LDA+U functional correctly describes the essential 
physics of the weakly correlated metals: (i) reducing the band dispersion by dressing of 
the one-particle excitation, and (ii) spin fluctuations near the quantum critical point. One 
functional, labeled FLL here, correctly describes the important physics in the limit of well 
localized electrons, and can be recommended in this case. The other functional, labeled 
AMF, is exact in a hypothetical material with the uniform orbital occupancies. Although 
neither functional accounts for the fluctuation effects, LDA+U may be useful, if applied 
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with a grain of salt, in moderately correlated metals. For this case, we propose a recipe that 
accounts for an incomplete localization and reduces to AMF or FLL in the appropriate limits. 
Finally, it is worth noting that in many correlated materials the spin-orbit interaction plays 
a key role. Since our a does not depend on spin, this prescription can be also formulated 
in terms of the full (41 + 2) x (4/ + 2) occupation matrix p and n = Tr(p)/(4Z + 2). Eq. 
(H) should be replaced with Tr(5p ■ 5p) = (4/ + 2)an(l — n). This formulation has another 
advantage in the case of a half-filled band, like in Gd, because in this limit it reduces to more 
physically meaningful in this case FLL, rather than to AMF as the nonrelativistic Eq. (U). 
This work is supported in part by ONR and by NSF (Grant No DMR-0071823). 
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